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Starting from the meson-meson Green function in 3+1 dimensional quenched lattice QCD we calculate poten- 
tials between heavy-light mesons for various light-quark mass parameters. For the valence quarks we employ the 
staggered scheme. The resulting potentials turn out to be short ranged and attractive. A comparison with a 
tadpole improved action for the gauge fields is presented. 
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1. INTRODUCTION 

For several decades nucleon-nucleon interac- 
tions have been parametrized by phenomeno- 
logical potentials. Substantial effort, employing 
purely hadronic degrees of freedom, has lead to 
meson-exchange potentials Attempts to take 
into account quark and gluon degrees of freedom 
with hybrid models have also been made @]. 

Today, QCD is believed to be the fundamen- 
tal theory of strong interactions. Thus it has be- 
come a challenge for theoretical nuclear physics 
to extract a nucleon-nucleon potential from first 
principles. In the low energy regime of QCD non- 
perturbative tools have to be used. This leads us 
to study the forces in systems of two hadrons on 
the lattice. 

Previous lattice calculations with static valence 
quarks have revealed an attractive force between 
two three-quark clusters (3). When dynamical 
quark propagators are used antisymmetrization 
and the exchange of valence quarks become pos- 
sible (J]. In 2+1 dimensional QED the potential 
between light mesons exhibits a repulsive hard 
core and is attractive at intermediate distances 
|U . An extension of this formalism to 3+1 dimen- 
sional QCD using a hopping parameter expansion 
for the quark propagators is reported in Q . 

In the present work we take another step to- 
ward the goal of calculating hadronic potentials 
from QCD. We study a system of two heavy-light 
mesons with dynamical quark propagators for the 



light valence quarks. Results from calculations 
with the Wilson action for the gauge field and a 
tadpole improved gauge field action are reported. 

2. THEORY 

2.1. Meson-meson correlator 

The one-meson field is a product of staggered 
Grassmann fields \ and % with a heavy and a 
light external flavor h and I, respectively, 

Mt) = xi(xt)xh($t) . (1) 

The meson- meson fields with relative distance r = 
y — x are then constructed by 



$ r -(t) = v- 1 z2 Sw-afeWM*) 



(2) 



The information about the dynamics of the 
meson-meson system is contained in the time cor- 
relation matrix 

CW'(Mo) = ($tW*r'(io))-(4W)(*?'(io)),(3) 

where ( ) denotes the gauge field configuration 
average. On the hadronic level C is a 2-point cor- 
relator of a composite local operator describing a 
molecule-like structure. Working out the contrac- 
tions between the Grassmann fields the following 
diagrammatic representation is obtained 



c = - = 



M 



(4) 



Each contribution to the correlator comprises 
the exchange of gluons. Thus, even diagram 
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leads to forces, whereas diagram cor- 
responds also to interaction mediated by the ex- 
change of the light valence quark. Denoting con- 
tractions of the Grassmann fields by 



X...X 



G,, 



with n = 1 . 



(5) 



we have for example 
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CM ~ (4,4,0^) = (GW(T®&Q(%%(6) 

with f= y — x, r' — y' — x' , and ~ stands for the 
sums and factors that carry over from (0) . Color 
indices are suppressed. The gauge configuration 
average is taken over the product of all four prop- 
agators G. The propagator of the light quark G^ 
is obtained from inverting the staggered fcrmion 
matrix with a random source estimator. A stan- 
dard conjugate gradient algorithm is used. The 
heavy-quark propagator is given by 



xt,xt a 
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where the phase factors T, M = (-l)^+^+ x ^/ a 
in the Kogut-Susskind formulation correspond to 
the Dirac matrices and k = (t — to) /a. In our 
calculations we set mh,a = 1. 

Since the heavy valence quarks are fixed in 
space the relative distance between the mesons 
is the same at the initial and final time of the 
propagation, r' = f. The effective ground-state 
energy of the meson-meson system W{r) can then 
be extracted from the large euclidean time behav- 
ior of C following quantum-mechanical reasoning 
for composite particles 0, 

C rr <t,t ) = Y,\^ n )\ 2 ^ En{t ' to) 

n 

~ c(r ) e -^(*0(*-*o) . (8) 
The residual meson-meson potential is 
V{r) = W{r)-2m, (9) 
with the mass 2m of two free mesons subtracted. 

2.2. Improved action 

Discretization errors due to finite lattice spac- 
ing a can be reduced by including terms of higher 



order in a with the action. It has been shown that 
at the classical level adding a term with six-link 
rectangular plaquettes U r t to the usual Wilson 
gauge field action removes 0(a 2 ) errors 0. A 
significant improvement is obtained by introduc- 
ing tadpole factors in the six- link term ||]. We 
use the improved action 



S[U] = PplJ2 \^eTr(l-U pl ) 



P i 



f3 r tY,o^eTr(l~U rt ), 



(10) 



where the first term is the Wilson action with 
four-link plaquettes U p i . The coupling parameter 



0T 



0PI 

20ul 



(11) 



and the mean link 
U o = (^ReTr(C/ p; )) 1/4 



(12) 



are determined self-consistently for a given fj. 
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3. RESULTS 



We considered a periodic 8 3 x 16 lattice. Each 
potential is the result of a measurement on 100 in- 
dependent gauge field configurations which were 
separated by 200 sweeps. The inversion of the 
fcrmion matrix was performed with 32 random 
sources. 

The effective energy W(r) was extracted from 
the correlator C by a four parameter (Levenberg- 
Marquardt) fit with the function 

C rr {t) = B(r) cosh [W{r){t -8a)] 

+ (-l)^ a B(r) cosh [W(r)(t - 8a)}. (13) 

The second term alternating in sign is a peculiar- 
ity of the staggered scheme. In order to improve 
the quality of the fits 15 time slices were used. 
The extracted potentials are thus between meson 
states which are contaminated by excitations. To 
be numerically consistent with W(r) at large dis- 
tances, the mass 2m of two noninteracting mesons 
was extracted from the square of the meson two- 
point function [C^] 2 . 
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Figure Oj: Potentials between heavy-light mesons obtained from simulations with conventional Wilson 
action compared to results with improved action. The potentials are attractive and become stronger with 
decreasing light-quark mass. Error bars are smaller when using tadpole improved action. Improvement 
was restricted to the gauge field action. 
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3.1. Wilson action 

For the update with the usual Wilson action we 
chose the inverse gluon coupling constant = 5.6. 
This corresponds to a lattice spacing of approxi- 
mately 0.2fm. The left column of Fig. |l| shows the 
resulting potentials for different values of the light 
valence quark mass mia = 0.05,0.1 and 0.2, re- 
spectively. In all cases the potential between two 
heavy-light mesons is short ranged and attractive. 
Calculations from inverse scattering theory pro- 
pose a similar shape for KK-potentials Q . Exper- 
imental scattering data contain inelastic channels 
which are not accounted for in our QCD calcu- 
lation. At distances r/a > 1 essentially no in- 
teraction energy can be resolved and the effective 
total energy approaches the value of two noninter- 
acting mesons 2m. With decreasing light-quark 
mass the interaction becomes stronger. 

3.2. Improved action 

Here the coupling was set to f3 p i — 7.6. This 
value corresponds to a w 0.2fm, as for the Wilson 
action ||. The right column of Fig. [l] shows the 
resulting potentials, again for different values of 
the light valence quark mass m;et = 0.05, 0.1 and 
0.2. Again we obtain attractive potentials which 
become deeper with decreasing quark masses. 
Those exhibit larger anisotropies than in the Wil- 
son case. The likely reason is that only the glu- 
onic action but not the staggered fermion matrix 
has been improved. In any case, it is interesting 
to observe that the shape of the potential remains 
stable. 

4. CONCLUSIONS 

We applied a practical method to extract effec- 
tive hadron-hadron potentials from lattice QCD 
to a system of two heavy-light mesons. For the 
light valence quark we employed dynamical quark 
propagators. The resulting potentials are short 
ranged and attractive. Smaller quark mass pa- 
rameters lead to stronger interaction. A prelim- 
inary analysis of simulations with improved ac- 
tion showed consistent results. Further insight 
into the interaction mechanism is expected from 
calculations with mesons consisting of two light 
valence quarks. 
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